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Abstract 

We propose a generalization of Grnnberg’s method of effective charges in 
which, starting with the effective charge for some dimensionless QCD observ¬ 
able dependent on the single energy scale Q, Ti{Q), we introdnce an infinite 
set of anxiliary effective charges, each one describing the snb-asymptotic Q- 
evolntion of the immediately preceding effective charge. The corresponding 
infinite set of conpled integrated effective charge beta-fnnction eqnations may 
be trnncated. The resnlting approximations for TZ{Q) are the convergents of 
a continued function. They are manifestly RS-invariant and converge to a 
limit equal to the Borel sum of the standard asymptotic perturbation series 
in asiO, with remaining ambiguities due to infra-red renormalons. There 
are close connections with Fade approximation. 
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QCD perturbation theory in its conventional formulation of a power series 
expansion in the renormalization group (RG) improved coupling, as(/i^), 
suffers from a number of defects. The series has a zero radius of convergence 
with coefficients exhibiting factorial growth, and consequently its sum can 
only be reconstructed by treating it as asymptotic to some function of the 
coupling. Customarily Borel summation is used for this reconstruction but 
of course this is not the only possible choice, and the method to be used is 
not uniquely specihed by the physical theory itself [|l|. Whilst ultra-violet 
(UV) renormalon singularities in the Borel plane cause no problems in the 
dehnition of the Borel sum, there are ambiguities resulting from infra-red (IR) 
renormalons which, it has been argued , are closely connected with, and can 
compensate, corresponding ambiguities in the operator product expansion 
(OPE) 0. 

In addition to these problems arising from the large-order behaviour of 
conventional perturbation theory, hxed-order perturbation theory is afflicted 
with the problem of renormalization scheme (RS) dependence. Thus hxed- 
order truncations of the series depend on the renormalization convention 
used to dehne as(/i^). A number of proposals for controlling or avoiding this 
RS-dependence problem have been advanced [3-6], but no consensus on the 
issue has been reached. 

In this paper we suggest a reformulation of QCD perturbation theory 
which avoids many of these problems. In particular it is convergent, notwith¬ 
standing remaining ambiguities due to IR renormalons, and is manifestly 
RS-invariant. It can be regarded as a generalization of Grunberg’s effective 
charge approach [^ and turns out to have a close relation to Bade approxi¬ 
mation 0, ID . 

We begin by considering a generic dimensionless QCD observable TZ{Q), 
which depends on the single energy scale Q (e.g. the R-ratio with Q 
the c.m. energy). By raising to a power and scaling we can always arrange 
that the formal perturbation series for TZ{Q) assumes the form, 

= a r\a^ -|- r2(i^ -|- ... - 1 - -|- .. . , ( 1 ) 

where a = Q!s(/i^)/7r is the RG-improved coupling . We shall refer to observ¬ 
ables dehned in this way as effective charges |Q, |[. 

Such effective charges satisfy several important properties. We shall con¬ 
sider SU(A^) QCD with Nf flavours of massless quark, and assume that the 
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first beta-function coefficient b = {llN—2Nf)/6 is positive. The first familiar 
property is Asymptotic Freedom (AF), which is equivalent to the statement 
that for any effective charge Ti{Q), 

lim 'Tl{Q) = 0 . (2) 

Q^oo 

The second property, which we shall refer to as “Asymptotic Scaling” (AS) , is 
perhaps less familiar, and will be central to the reformulation of perturbation 
theory we are proposing. The hrst step is to dehne a universal QCD scaling 
function. 


T{x) = + l/cxf^ , 


where c is the second (universal) beta-function coefficient. 


(3) 


5 3 ■ 


(4) 


with Ca = N, Cf = {N'^ — 1)/2A^ , SU(A^) Casimirs. Then AS corresponds 
to the statement that for any effective charge 77(Q), 


hm QTinm = Kn , (5) 

Q^oo 

where Kn is a scaling constant with dimensions of energy which depends 
on the observable. Given sufficiently large values of Q this property can 
evidently serve as a test of QCD, but since K-ji is not universal it cannot 
usefully be applied at hxed values of Q . However if the next-to-leading order 
(NLO) perturbative coefficient r = rj^‘^(/r = Q) has been calculated (in the 
MS scheme with scale /i^ = for instance), then K-ji can be converted 
into a universal scaling constant via the exact Celmaster-Gonsalves [Q 
relation, 

At^c ^ = Ajjg . (6) 

Use of different subtraction procedures results in different constants which, 
however, are still universal. We stress that r is independent of Q. Thus we 
arrive at the property of “Universal Asymptotic Scaling” (UAS), 

lim 0.;r(75(Q))e--/i' = A (7) 

Q—>00 
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This property can now be nsed to test QCD at fixed Q by looking at the 
scatter in for varions observables 0. At finite Q UAS will 

be violated and one will have 


Qr{n{Q))e-->'’ = Ajjj(l + r,K(Q)) , (8) 

where ri'ji{Q) represents the non-scaling snb-asymptotic effects which vanish 
as Q —> cxo . The crncial observation is that r]Tz{Q) has the formal pertnrba- 
tion series 


rmiQ) = j(a + 0(a")) , (9) 

where p 2 is the next-NLO (NNLO) effective charge (EC) beta-fnnction coef¬ 
ficient 1^, 1^, given by the RS-invariant combination , 

p2 = C2 + r2- riC - ri^ , (10) 

with C 2 the three-loop beta-fnnction coefficient. Thns pn is proportional to 
another effective charge, which will itself satisfy the UAS property, and will 
have its own non-scaling contribntion proportional to yet another effective 
charge, etc. This self-similar constrnction based on the exact property of 
UAS provides the reformnlation of QCD pertnrbation theory which we are 
proposing. 

The constrnction can nsefnlly be regarded as a generalization of Grnn- 
berg’s effective charge approach |^, and it will help to illnminate the 
discnssion so far if we briefly review that formalism. We refer the reader to 
reference [6] for fnll details. 

The key ingredient in the EC approach is the EC beta-fnnction p{lZ{Q)) 
which determines the Q-evolntion of the observable TZ{Q), 


dn{Q) 

d\nQ 


-hpinm 

-b{n^ + c7^3 + + ... + + ...), (11) 


where the pn are Q-independent RS-invariant combinations of the r,, q {i < 
n) (see eq.(lO) for P 2 ). p{lZ{Q)) can be regarded as an observable to be 
reconstrncted from the measnred rnnning of 7^(Q) |P]. Integrating np eq.(ll) 
and imposing Asymptotic Freedom as a bonndary condition yields. 
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F(7^(Q)) 


= bln^-T^^^dx 

An Jo 

= b\n-^ - ApoiQ) , 

J-^n 


1 1 
p{x) x‘^{l + cx) 


(12) 


where 

F{x) =—h\n [x) = — + chi {—^^ , (13) 

X \1 + cxj 

and with At^ a constant of integration. Assnming AF the AS property in 
eq.(5) then follows directly on rearranging eq.(12) and taking the Q ^ oo 
limit, in which Apo(Q) —> 0 . Further , as demonstrated in ref.[6], one can 
then identify hhi{Q/A ti) with the NLO RS-invariant po(Q) Bi. 

61nA = p„(Q) = Mn^-r, (14) 

A7^ J-y-Ms 

from which the Celmaster-Gonsalves relation in eq.(6) follows. The sub- 
asymptotic pniQ) contribution in eq.(8) is equivalent to (exp(Apo/&) — 1), 
with Apo(Q) as dehned in eq.(12). 

In this notation the violation of UAS is controlled by Apo(Q) • Expanding 
the integrand in eq.(12) as a power series in x and integrating term-by-term 
we obtain 


'^Po(Q) — P2’^+(P3~2p2c)^--|-(p4 —2cp3-|-2c^P2 + 2c^—+ - • • • (15) 

We can then write Apo(Q) = P 2 TiJ^\Q)-, where 71A\Q) is another effective 
charge whose perturbative expansion follows on substituting TZ = a + rioJ -|- 
.. . in eq.(15), 


'RA\Q) = a + {ri + — c)o? + , (16) 

2p2 

from which the coefficients ..., can be obtained. TZA\Q) will satisfy 

eq.(12) with po(Q)) 2kpo(Q) replaced by p^q\Q) and Apg^^((5) where 
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Po\Q) = b\n-^-r\^^^‘‘^ = b\n-^-r - ^ + c. ( 17 ) 

Writing Ap^q\Q) = with VP\Q) yet another effective charge 

the self-similar construction continues. 

Thus perturbation theory is reformulated as an inhnite set of effective 
charges (Q), .. .,77*^"'^((5),. • satisfying the set of coupled equa¬ 

tions, 


F(nQ)) = Po{Q)-P2'li''\Q) 
F(K<‘'(0)) = p|i‘’W)-p5‘’K®(0) 
f(R<">(Q)) = p?'(Q)-pfKW(Q) 

F(n'"HQ)) = py(Q)-p^"’K'’‘+‘'(Q) 


(18) 


Po"* involves bln(Q/r, and the pi; EC RS invariants for k < 2n + 1 (we 
dehne pi = c). Thus it is determined given a complete perturbative 

calculation for Tl{Q) . p'^^ involves the pk ior k <2n + 2 and so is determined 
given a calculation for TZ{Q). 

The effective charge describes the violation of UAS of 77((5), 

TZZ^\Q) describes the scaling violation of TZ^^^Q) etc. If we assume that 
TZ^'^\Q) scales and set lZb^^^\Q) = 0, then the truncated set of equations in 
eqs.(18) can be solved for 7Z{Q), given ln((5/Ajg;5), the NLO MS coefficient 
r, and the Pq \ p 2 \ fhat is given a complete N^”+^LO perturbative calculation 
for TZ{Q). 

Equivalently we can iteratively eliminate TZ^^\ TZZ^\ ... from eqs.(18). If 
we dehne G{x) to be the inverse function of F{x) in eq.(13), that is G{F{x)) = 
x, then we obtain the continued function representation for 7Z{Q), 


K(Q) = G(p„ - p^GipP - p^‘'G(pf - pf G(...)))...) . (19) 

Equation (19) is our main result. The hrst NLO approximation G{po) 
is precisely the conventional NLO perturbative approximation in the EC 
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scheme. Subsequent approximations differ from conventional fixed-order per¬ 
turbation theory. An interesting feature of the approach is that each extra 
level in the construction of eqs.(18) adds two further orders in perturba¬ 
tion theory. One can interpolate by using G{po — p 2 G{po)) as the NNLO 
approximation, and then the N^LO approximation is G{po — P 2 G{p^q'^)), 
G{po — P 2 G{pI^'^ — interpolates the N^LO, etc. 

Our claim is that if only ultra-violet (UV) renormalon singularities are 
present then this sequence of approximations converges, and the limit is 
precisely the Borel sum of the divergent asymptotic conventional perturba¬ 
tion series in as{p^)- If infra-red (IR) renormalons are present in addition 
then strict convergence is prevented [§]. Since by construction eqs.(18) are 
a set of exact equations this indicates that these equations are insufficient 
to determine TZ{Q) and correspondingly extra information must be added, 
presumably in the form of further OPE terms. 

We shall now motivate this convergence claim by considering the simpli- 
hed case of QCD with c = 0. In this limit F{x) = G{x) = 1/x and the 
continued function of eq.(19) becomes a continued fraction. 


nQ)^ 


P2 




( 20 ) 


IPo Ipo Ipo 

As we shall now show the successive convergents of this continued fraction are 
the diagonal [n/n] Fade approximants of the original perturbation series for 
7?.(Q) in eq.(l), in the RS where C 2 = C 3 = ... = = ... = 0, the so-called 

‘t Hooft scheme [|^. This is in accord with the recent observation by Gardi 
that in the “large- /3o” limit diagonal Fade approximants are RS-invariant 


The diagonal [n/n] Fade approximants are dehned by 


7^[n/n] 


a ^ 2 ®^ -l- ... -|- 
1 -|- B\a -j- ... -|- 

a rlO? -\- r 2 a/ -1- ... -|- T 2 n—io/^^ + , ( 21 ) 


where the Ai and Bi coefficients are uniquely hxed by demanding that the 
rj,i < 2 n — 1 are reproduced on expanding the quotient. 

Corresponding to the original perturbation series in eq.(l) one can dehne 
the standard Stieltjes continued fraction for 71{Q), 
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( 22 ) 


K«3) = + ^>1 + 


1 


1 


where Ki = 1, K 2 = —ri, = ri — r 2 /ri,.... The successive convergents of 
eq.(22) are respectively the diagonal [n/n] and off-diagonal [n/n + 1] Fade 
approximants to 1Z{Q). We can also dehne the so-called “associated” 
Jacobi form of the continued fraction 


or 
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7^(Q) = 


Kia\ 


K2K3a^ 




(23) 


1(1 J- K2C1) |1 J- (-K^s J- -^"4)0 |1 -|- {K^ Koja 

The convergents of eq.(23) are the diagonal [n/n] Fade approximants of TZ. 
If we scale the partial quotients in eq.(23) by a we obtain 


7^(Q) = 


Ki\ 


KoK, 


K.K. 


|(l/a)-|-iF2 |(l/a)-|-iFs-|-7^4 \{l/a) + K^ + KQ 


• • (24) 


In the ‘t Hooft scheme (c 2 = C 3 = ... = 0), and with c = 0 , = [1/a)—r^i\ 

and so on comparing eq.(24) with eq.(20) we can identify = (1/a) -|- 
J^ 2 n+i + J^ 2 n +2 and p^”^ = iF 2 n-i -2 J^ 2 n+ 3 , and we see that indeed the successive 
(manifestly RS-invariant) convergents of eq.(20) are precisely the successive 
convergents of the Jacobi form of the continued fraction- the diagonal 7 ^[n/? 7 ,] 
Fade approximants. 

The convergence of the diagonal Fade approximants is assured if the 
perturbation series for 1Z{Q) is asymptotic to a “Stieltjes function” OH. 
and if this is the case this function is identical to the Borel sum of the series, 
if it exists . 


A Stieltjes function is of the form 


fC 

71(a) = a 

Jo 


dt- 


W(t) 


(25) 


(1 -|- ta) ’ 

where W(t) is nonnegative throughout the range of integration. For such a 
function the continued fraction coefficients Kn are all nonnegative 
A relevant example, motivated by UV renormalon singularities 
consider a simple pole at 2 ; = —1 in the Borel plane. 
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is to 
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( 26 ) 


roc 

71(a) = dzi 
Jo 


-zla_ 




= a 


dt- 


(1 + z) Jo (1 + ta) 

so that the non-negative Stieltjes weight function in eq.(25) is W(t) = e“*. 
This corresponds to perturbative coefficients = (—l)”n! which results in 
the (nonnegative) continued fraction coefficients = 1, K 2 n = K 2 n+i = 
n,n > 1, from which we can obtain = (1/a) -|- 2n -|- 1 and = {n + 1)^ 
for the coefficients in eqs.(18). The continued fraction in eq.(20) then con¬ 
verges to the well-defined Borel sum in eq.(26). Branch point UV renormalon 
singularities can also be reduced to Stieltjes form by a change of variables. 
A subtlety is that if one changes the RS , and hence the definition of a, the 
series may no longer be Stieltjes, equivalently the Kn are RS-dependent and 
will not all be nonnegative in a general RS. To prove convergence, however, 
it is only necessary that the series is Stieltjes in one particular RS, since the 
convergents are RS-invariant. 

By analyzing the large-n behaviour of eqs.(18) using 2n and p^"^ ~ 

v? one can infer that for the above example 7lJ'^\Q) ~ 1/n, so as n —>■ 
CX), 7&"\Q) —0, and successive violations of UAS become ever smaller, 
underwriting the convergence of the truncations of eqs.(18). 

To conclude, we regard the important feature of this approach to be its 
physical motivation based on iterating an exact property of massless QCD, 
which we have termed Universal Asymptotic Scaling, and which is dehned in 
eq.(7). Thus if the inhnite set of coupled equations in eqs.(18) has a solution 
this determines TZ{Q) ■ It is incidental that this approach can reproduce 
the result of applying Borel summation to the original divergent asymptotic 
series. For instance if = (—1)"^(2?7,)! the Borel method is inapplicable, 
nonetheless the set of coupled equations eqs.(18) in this case generates a 
continued fraction (assuming c = 0) which converges to a unique Stieltjes 
function . 

We plan to give a more complete discussion of the convergence of the 
approach, and the manner in which ambiguities induced by the presence of 
IR renormalons manifest themselves in this language, in a future work. The 
links with Bade methods are also interesting and should be explored further. 
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